Critical Behavior in Light Nuclear Systems: Experimental Aspects 
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An extensive experimental survey of the features of the disassembly of a small quasi-projectile 
system with A ~ 36, produced in the reactions of 47 MeV/nucleon 40 Ar + 27 Al, 48 Ti and 58 Ni, 
has been carried out. Nuclei in the excitation energy range of 1-9 MeV/u have been investigated 
employing a new method to reconstruct the quasi-projectile source. At an excitation energy ~ 5.6 
MeV/nucleon many observables indicate the presence of maximal fluctuations in the de-excitation 
processes. These include the normalized second moments of the Campi plot and normalized variances 
of the distributions of order parameters such as the atomic number of the heaviest fragment Z max 
and the total kinetic energy. The evolution of the correlation of the atomic number of the heaviest 
fragment with that of the second heaviest fragment and a bimodality test are also consistent with 
a transition in the same excitation energy region. The related phase separation parameter, S p , 
shows a significant change of slope at the same excitation energy. In the same region a A-scaling 
analysis for of the heaviest fragments exhibits a transition to A = 1 scaling which is predicted to 
characterize a disordered phase. The fragment topological structure shows that the rank sorted 
fragments obey Zipf's law at the point of largest fluctuations providing another indication of a 
liquid gas phase transition. The Fisher droplet model critical exponent r ~ 2.3 obtained from the 
charge distribution at the same excitation energy is close to the critical exponent of the liquid gas 
phase transition universality class. The caloric curve for this system shows a monotonic increase 
of temperature with excitation energy and no apparent plateau. The temperature at the point of 
maximal fluctuations is 8.3 ± 0.5 MeV. Taking this temperature as the critical temperature and 
employing the caloric curve information we have extracted the critical exponents /3, 7 and a from 
the data. Their values are also consistent with the values of the universality class of the liquid 
gas phase transition. Taken together, this body of evidence strongly suggests a phase change in an 
equilibrated mesoscopic system at, or extremely close to, the critical point. 



PACS numbers: 25.70.Pq, 24.60.Ky, 05.70.Jk 
I. INTRODUCTION 



Probing the liquid gas phase transition of finite nuclei 
is an important topic in nuclear physics since it should 
allow investigation of the nuclear equation of state and 
clarify the mechanism by which the nucleus disassem- 
bles when heated. This phase transition is expected to 
occur as the nucleus is heated to a moderate tempera- 
ture so that it breaks up on a short time scale into light 
particles and intermediate mass fragments (IMF). Most 
efforts to determine the critical point for the expected liq- 
uid gas phase transition in finite nuclconic matter have 
focused on examinations of the temperature and excita- 
tion energy region 0, Q where maximal fluctuations in 
the disassembly of highly excited nuclei are observed [3j . 
A variety of signatures have been employed in the identifi- 



cation of this region 0, IE IE 00 an d several publications 
[EUEm] have reported the observation of apparent crit- 
ical behavior. Fisher Droplet Model analysis have been 
applied to extract critical parameters |12(. The derived 
parameters are very close to those observed for liquid-gas 
phase transitions in macroscopic systems 13] . Data from 
the EOS [3 and ISiS HE Eg collaborations have been 
employed to construct a co-existence curve for nucleonic 
matter |l2j . Interestingly, the excitation energy at which 
the apparent critical behavior is seen is closely correlated 
with the entry into the plateau region in the associated 
caloric curve [lj. Although implicit in the Fisher scal- 
ing analyses is the assumption that the point of maximal 
fluctuations is the critical point of the system 0, ^} , 
other theoretical and experimental information suggest 
that the disassembly may occur well away from the criti- 
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cal point Q, H E EE E] . In addition, recent lattice 
gas calculations indicate that the Fisher scaling may be 
observed at many different densities [22 > raising doubts 
about previous critical point determinations. Further, 
applications of A-scaling analysis indicate that the obser- 
vation of power-law mass distributions jU |3| , although 
necessary, is not sufficient to identify the true critical 
point of the system. We note also that while the role of 
the long range Coulomb interaction in determining the 
transition point has received considerable theoretical and 
experimental attention E E E E E E a number 
of questions remain as to the appropriate way to deal 
with the complications it introduces. 

In this paper we report results of an extensive investi- 
gation of nuclear disassembly in nuclei of A ~ 36, excited 
to excitation energies as high as 9 MeV/nucleon. To our 
knowledge, this is the smallest system for which such an 
extensive analysis has been attempted. An earlier brief 
report on some aspects of this work appeared in Physical 
Review Ej. While investigating a smaller system takes 
us farther from the thermodynamic limit, several theoret- 
ical studies indicate that phase transition signals should 
still be observable E,E] in small systems. The choice of 
a lighter system for investigation has the advantage of re- 
ducing the Coulomb contributions. Earlier investigations 
in this mass region have provided valuable insi ghts into 
the binary reaction mechanism |33l 133. Ei 1361 1371 |38| , 
multifragment emission as a function of excitation en- 
ergy Ha E E E E B E], the emission time-scale 
and emission sequence of light particles E> E> E| an d 
collective flow behavior 0| etc. 

Applying a wide range of methods we find that the 
maximum fluctuations in the disassembly of A ~ 36 oc- 
cur at an excitation energy of 5.6±0.5 MeV and a tem- 
perature of 8.3±0.5 MeV. At this same point, the crit- 
ical exponents describing the fragment distributions are 
found to be very close to those of the universality class 
of the liquid gas phase transition. 

These observations do not guarantee critical behavior 
has been reached, however, in contrast to experimental 
results for heavier systems QJ we also find that the caloric 
curve for A ~ 36 does not exhibit a plateau at the point 
of maximum fluctuations. Taken together, the observa- 
tions strongly suggest a phase change in an equilibrated 
mesoscopic system at, or extremely close to, the critical 
point. 

The paper was organized as follows: in Sec. II we de- 
scribe the set-up of our experiment; Section III presents 
a new method for reconstruction of the quasi-projectile 
source, QP; Section IV discusses some general features of 
the reconstructed QP; Section V explores the evidence for 
critical behavior in the disassembly of the QP. In Section 
VI, we discuss the caloric curve of the QP In Section VII 
we use the scaling theory to derive the critical exponents 
of the transition. All those values are found to be con- 
sistent with the universality class of the liquid gas phase 
transition. Conclusions are presented in Section VIII. 



TABLE I: NIMROD Charge Particle Array 
Ring Angle (deg) No. of Segments Solid Angle(src) 
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II. EXPERIMENTAL SET-UP AND EVENT 
SELECTION 

Using the TAMU NIMROD (Neutron Ion Multide- 
tector for Reaction Oriented Dynamics) E] and beams 
from the TAMU K500 super-conducting cyclotron, wc 
have probed the properties of excited projectile-like frag- 
ments produced in the reactions of 47 MeV/nucleon 40 Ar 
+ 27 A1, 48 Ti and 58 Ni. Earlier work on the reaction 
mechanisms of near symmetric collisions of nuclei in the 
20 < A < 64 mass region at energies near the Fermi 
energy have demonstrated the essential binary nature of 
such collisions, even at relatively small impact parame- 
ters As a result, they prove to be very useful in 
preparing highly excited light nuclei with kinematic prop- 
erties which greatly simplify the detection and identifi- 
cation of the products of their subsequent de-excitation 
0. 

The charged particle detector array of NIMROD, 
which is set inside a neutron ball, includes 166 individual 
Csl detectors arranged in 12 rings in polar angles from 
~ 3° to ~ 170°. Eight forward rings have the same ge- 
ometrical d esig n as the INDRA detector, but have less 
granularity 52]. The angles, number of segments in each 
ring and solid angle of each Csl segment are given in 
Table I. 

In these experiments each forward ring included two 
super-telescopes (composed of two Si-Si-Csl detectors) 
and three Si-Csl telescopes to identify intermediate mass 
fragments. The Csl detectors are Tl doped crystals read 
by photo-multiplier tubes. A pulse shape discrimination 
method using different responses of fast and slow compo- 
nents of the light output of the Csl crystals is employed 
to identify particles 53] . In the Csl detectors Hydrogen 
and Helium isotopes were clearly identified and Li frag- 
ments are also isolated from the heavier fragments. In the 
super-telescopes, all isotopes with atomic number Z < 8 
were clearly identified and in all telescopes particles were 
identified in atomic number. 
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The NIMROD neutron ball, which surrounds the 
charged particle array, was used to determine the neutron 
multiplicities for selected events. The neutron ball con- 
sists of two hemispherical end caps and a central cylin- 
drical section. The hemispheres arc upstream and down- 
stream of the charged particle array. They are 150 cm in 
diameter with beam pipe holes in the center. The cen- 
tral cylindrical sections 1.25m long with an inner hole 
of 60 cm diameter and 150 cm outer diameter. It is di- 
vided into four segments in the azimuthal angle direction. 
Between the hemispheres and the central section, there 
are 20 cm air gaps for cables and a duct for a pumping 
station. The neutron ball is filled with a pseudocumene 
based liquid scintillator mixed with 0.3 weight percent 
of Gd salt (Gd 2-ethyl hexanoate). Scintillation from a 
thermal neutron captured by Gd is detected by five 5-in 
phototubes in each hemisphere and three phototubes in 
each segment of the central section. 

The correlation of the charged particle multiplicity 
(M cp ) and the neutron multiplicity (M n ) was used to 
sort event violence. In Fig. 2] lines indicate the event 
windows which have been explored. Roughly speaking , 
the more violent collisions correspond to those with the 
highest combined neutron multiplicity and charged par- 
ticle multiplicity (M cp )(Binl). This can be seen in the 
excitation energy distribution of the QP in Fig. (the 
determination of excitation energy will be explained in 
the following section). For Binl and Bin2, the average 
excitation is ~ 4 MeV/nucleon and the E*/A distribu- 
tion extends to 9 MeV. Since the goal of the present work 
was to explore the disassembly of highly excited QP, we 
have used the data from Binl and Bin2 together in the 
present work. For the more peripheral bins, i.e. Bin4 and 
Binb, however, there is apparent event mixing in the up- 
per range of excitation energy. In that case, events com- 
ing from a particular excitation energy can be distributed 
over several experimentally reconstructed excitation en- 
ergy bins. This judgment is supported by the data plot- 
ted in Figure 3, where we plot, for the QP formed in 
40 Ar + 58 Ni reactions, the total multiplicity of charged 
particles (Fig. |3Ja)). the QP normalized charge number 
of heaviest fragment (Fig. E^b)), the effective Fisher's 
power-law parameter r e ff (Fig. EJc)) and Zipf's law pa- 
rameter £ (Fig. Eld)) as a function of the excitation en- 
ergy in different centrality bins (for detailed explanations 
of the physical quantities plotted, see the following sec- 
tions). As seen in the figure, Binl and Bin2 display es- 
sentially identical behavior and the values do not depend 
the selected centrality bin. However, data for the periph- 
eral bins (Binb and BinA) deviate significantly from the 
data for Binl and Bin2 in the upper range of the re- 
constructed excitation energy. These deviations indicate 
event mixing in peripheral collisions and raises questions 
about the validity of the excitation energy determination 
in the upper excitation energy range for peripheral col- 
lisions. For the intermediate Bin3, results are close to 
those of Binl and Bin2. To minimize the effects of pos- 
sible event mixing and realize our goal of exploring the 
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FIG. 1: (Color online) Two dimension plot for Mcp vs Mn as a 
selector of collision centrality in the 40 Ar + 58 Ni reaction. 
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FIG. 2: (Color online) The distribution of excitation energy for 
different centrality bins for the QP formed in 40 Ar + 58 Ni. 



disassembly of highly excited QP, we choose only data 
for Binl and Bin2 and combine those results for further 
analysis. 

Given the limitation of IMF identification in NIMROD 
detectors which do not have Si associated with them, 
the events with complete, or near complete QP detection 
have to be isolated before the analysis of the QP fea- 
tures can proceed. In the following section we describe 
the techniques of QP reconstruction and event selection 
which we have employed . 



III. A NEW METHOD OF 
QUASI-PROJECTILE RECONSTRUCTION 

Intermediate energy heavy ion collisions are compli- 
cated processes in which the roles of the mean field 
and nucleon-nucleon interactions may both be important. 
Many reactions manifest the mixed features of both the 
low energy deep inelastic scattering mechanism and a 



4 



8 
7 
6 

. 5 
S 4 

3 

2 

10 



3= 





8 2 4 

<E*/A> (MeV) 



FIG. 3: (Color online) Comparison of some physical quantities as a 
function of excitation energy bin in different centrality bins for the 
QP formed in 40 Ar + 58 Ni. (a) The total multiplicity of charged 
particles; (b) The QP normalized charge number of the heaviest 
fragment; (c) The effective Fisher's power-law parameter and (d) 
Zipf's law parameter. The deviations in higher excitation energy 
of Bin5 and BinA from Binl and Bin2 can be attributed to the 
event mixture in the peripheral collisions. For detailed explanation 
of the physical quantities, see the following sections. 



high energy participant-spectator mechanism. 

It is well known that the laboratory frame kinetic en- 
ergy spectra of most light ejectiles can be reproduced 
with the assumption of emission from three different 
sources: a Quasi- Projectile (QP) source, an intermediate 
velocity or Nucleon-Nucleon (NN) source and a Quasi- 
Target (QT) source. To better understand origins of the 
emitted particles the ideal situation would be to have 
the ability to attribute each particle to its source on an 
event-by-event basis. However the spectral distributions 
from the different sources overlap significantly, making 
such an attribution not possible. Previous techniques to 
reconstruct QP have included identifying high velocity 
components of the QP l55l | or treating only particles 
emitted in the forward hemisphere in the projectile frame 
and then assuming identical properties for particles emit- 
ted in the backward hemisphere in order to recreate the 
QP source |33|, |3J|- Such a technique is limited in its 
application and not suited to situations in which fluctu- 
ations are to be investigated. 

For this work we have developed a new method for 
the assignment of each light charged particle (LCP) to 
an emission source. This is done with a combination of 
three source fits and Monte-Carlo sampling techniques. 
We first obtain the laboratory energy spectra for different 
LCP at different laboratory angles and reproduce them 



using the three source fits. In the laboratory frame, the 
energy spectra of LCP can be modelled as the overlap 
of emission from three independent moving equilibrated 
sources, i.e. the QP, NN and QT sources. For evaporation 
from the QT source, we take |56( 
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where Ei ab and T s are respectively the laboratory energy 
and apparent slope temperature. Mi are multiplicities. 
In the above formula, the Coulomb barrier is considered 
to be in the QT source frame, in this case, E' and E" 
are defined as 
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where v s is the magnitude of the source velocity and is 
taken along the beam direction. 9 is the angle between 
the source direction and that of the detected LCP. 

For the LCP from QP and NN, we take the Coulomb 
barrier in the laboratory frame [57j ■ For QP, we assume 
the surface emission form 
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The total energy distribution is the sum over the QP, 
QT and NN component. 

Fig. 0] shows examples of the three source fits for 
deuterons and tritons in the second most violent bin 
(Bin2). From these fits we know the relative contribu- 
tions from the of QP, NN and QT sources. Employing 
this information to determine the energy and angular de- 
pendent probabilities we analyze the experimental events 
once again and, on an event by event basis, use a Monte 
Carlo sampling method to assign each LCP to one of the 
sources QP, or NN, or QT. For example, the probability 
that a certain LCP i (i.e. p, d and t etc.) will be assigned 
to the QP source is 



Prob® p (Ei ab ,9, 
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FIG. 4: The three source fits for deutrons (left panels) and tritons 
(right panels) in (Bin2) events of 40 Ar + 58 Ni reaction. The lab 
angle is 6.4°, 18.2°, 32.1°, 64.2° and 120°, respectively, from top 
to bottom. The meanings of the symbols and lines are depicted in 
the bottom right panel. 



To illustrate the results of such a procedure, we show, 
in Fig. [SJ the velocity contour plots for protons to 
Lithium associated with the highest multiplicity windows 
in the 40 Ar + 58 Ni reaction. 

The panels on the left represent the data before any 
selection, i.e. the raw data with contributions from all 
emission sources. Obviously the particles are of mixed 
origin and it is difficult to make a meaningful physical 
analysis. Since we are interested in the QP source, we 
show, in the middle panels, the velocity contour plots 
for particles assigned to the QP source using the above 
reconstruction method. As expected from the technique 
employed, the results exhibit clean, nearly spherical, ve- 
locity contours, corresponding to isotropic emission in 
the rest frame of QP source. 

The projected parallel velocity distributions are de- 
picted in the right panels of Fig. El The solid histograms 
represent the total distribution and the histogram with 
the hatched area represents the contribution from the 
QP source. The peak velocity of this QP contribution 
is close to the initial projectile velocity although some 
dissipation is evident. 

Intermediate mass fragments, IMF, with Z>4 were 



FIG. 5: (Color online) The velocity contour plots for the light 
charged particles in violent events (Bin2) of 40 Ar + 58 Ni reac- 
tion. From top to bottom, protons, deutrons, tritons, 3He, a and 
Lithium ; the left panel shows the total velocity contour plot, the 
middle column depicts the velocity contour plots of particles from 
the QP source and the right column presents the corresponding dis- 
tributions of the parallel velocity from the total contribution (solid 
histogram) and from the QP source (hatched area). See details in 
text. 



identified in the telescope modules of NIMROD. For such 
ejectiles we have we have not used such fitting techniques. 
Rather we have used a rapidity cut (> 0.65 beam rapid- 
ity) to assign IMF to the QP source. We also checked the 
sensitivity of the above rapidity cut to the results, eg., us- 
ing > 0.55 or > 0.75 beam rapidity, there are only minor 
changes for source mass, excitation energy and tempera- 
tures etc, within ~ 10% error bars. This has no influence 
on any conclusions we draw in this article. Of course, 
this is expectable for such a binary-dominated reaction 
mechanism. Once we have identified all LCPs and IMFs 
which are assumed to come from the QP source, we can 
reconstruct the whole QP source on an event-by-event 
basis. 



6 




FIG. 6: (Color online) The correlation of total charge number of 
QP (Zqp) and the total multiplicity of charged particle (Mqp) for 
the violent events of 40 Ar + 58 Ni reaction. The region above the 
dotted line is used to reconstructed the quasi-projectile. Note that 
the scale of z-axis is logarithmic. 

Fig. shows the two dimensional plot of total charged 
particle multiplicity (Mqp) and total atomic number 
{Zqp) for the reconstructed QP source. Due to the lim- 
ited geometrical coverage of the telescopes, the efficiency 
of detection for nearly complete QP events is low. Note 
that the scale of z-axis of the figure is logarithmic. 

To select the nearly complete QP events, we choose 
events with Zqp > 12 (i.e. as good events. The part of 
the distribution above the line in Fig. 0) corresponds to 
the accepted region of violent collisions for the 40 Ar + 
58 Ni reaction. The reconstructed good events in that re- 
gion account for 4.3% of the total events for the violence 
bins selected (Binl and Bin2). For the reactions 40 Ar 
+ 48 Ti and 40 Ar + 27 A1, a similar portion of the good 
central events has been collected to make the same anal- 
ysis. Totally, 28000, 54000 and 56000 good QP events 
have been accumulated to make the following analysis 
for 40 Ar + 58 Ni, 48 Ti and 27 Al reactions, respectively. 
For this analysis the velocity of the QP source was de- 
termined, event by event, from the momenta of all QP 
particles. 



IV. GENERAL PROPERTIES OF THE 
EXCITED QP 

After the reconstruction of the QP particle source, the 
excitation energy was deduced event-by-event using the 
energy balance equation [H^ , where the kinetic energy of 
charged particles (CP), the mass excesses and the (unde- 
tected) neutron contributions were considered, i.e., 

M tot 

E* =Y, E i in ( CP ) + 2 M « T + Q ( 9 ) 

i=l 

where Q, the mass excess of the QP system is determined 
from the mass difference between the final QP mass, Aqp 



and the sum over the masses of the detected particles of 
the reconstructed QP, Y*ti M CP )- MQ P ) is deter- 
mined from the total reconstructed charge of the QP, as- 
suming the QP has the same N/Z as the initial projectile 
and Ai(CP) is the mass of each detected charged particle, 
which was calculated from the measured Zi(CP) thr ough 
the numerical inversion of EPAX parameterization |59| , 
except for Z = 1 and 2 for which experimental mass iden- 
tification was achieved. The neutron multiplicity M n 
was obtained as the difference between the mass number 
(Aqp) of the QP and the sum of nucleons bound in the 
detected charged particle, i.e., M n = Aqp — Y^Ai(CP). 
E^ m (CP) is the kinetic energy of the charged particles 
in the rest frame of QP. The contribution of the neutron 
kinetic energy was taken as 3/2M„T with an assumed T 
that is equal to that of the protons. As our detector is not 
100% efficient wc corrected observed events (on the av- 
erage) for undetected mass and energy. For a particular 
excitation energy bin the missing multiplicity for a given 
ejectile is the difference between the multiplicity derived 
from the source fit and the average detected ejectile mul- 
tiplicity for events in the acceptance window. Assuming 
that missed particles have the same average kinematic 
properties as the detected particles of the same species 
allows the appropriate corrections to be made. Since al- 
most complete projectile-like species were selected ini- 
tially the missing particles were usually protons. Using 
these techniques we find that the average QP has a mass 
of 36 and a charge of 16. 

Assuming the mean mass of missing particles in a given 
E* / A window is equal to A A , the contribution of missing 
excitation energy AE* can be approximated as AE* = 
A A- yj X°p P ' wnere ^meas 1S the excitation energy before 
the correction and ^2 Alcp is the sum of the masses of 
LCP ( A< 7) and neutrons in the same E* /A window. 
Thus the real excitation energy should be £J, eos + AE* . 
Filtering results of AMD-GEMINI calculations [13 by 
applying experimental acceptances leads to very similar 
corrections to those employed. 

Fig. [7] depicts normalized excitation energy distribu- 
tions for the selected QP events in 40 Ar + 27 Al (open 
circles), 48 Ti (open triangles) and 58 Ni (solid squares) 
for Binl + Bin2. These distributions are very similar. 
For violent collisions, the highest excitation energy of the 
QP can reach 9 MeV/nucleon. In the following analysis, 
we will generally separate the excitation energy distribu- 
tions into 9 windows, as shown by the slices in Fig. [3 For 
simplicity, we sometimes identify these E* /A windows as 
Excl through Exc9. 

Fig. [S] shows the total multiplicity distribution of 
charged particles in 9 excitation energy windows. For 
the quasi-projectiles formed in Ar induced reaction with 
different targets, the distributions keep the nearly same 
which is a reasonable results thanks of a clean QP recon- 
struction technique. Fig. [5] shows average multiplicity of 
LCP as a function of excitation energy. For p, d, t and 
3 He, the multiplicity rises monotonically but for a and 
Li, the multiplicities peak at E* /A near 6 MeV/nucleon 



7 



10" 



c 

Q) 9 

>10 

T3 



»10 3 



10 



I 



1 23456789 10 
EVA (MeV) 



FIG. 7: (Color online) Excitation energy distribution of QP formed 
in 40 Ar + 27 Al (open circles), 48 Ti (open triangles) and 5s Ni (solid 
squares). Dotted lines indicate the selected excitation energy bins. 
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FIG. 8: (Color online) The total multiplicity distribution of 
charged particles from the QP systems formed in 40 Ar + 27 Al 
(open circles), 40 Ti (open triangles) and 58 Ni (solid squares). 



and then drop at higher excitation energy. This behav- 
ior is similar to the rise and fall behavior of IMF yield 
observed in many previous multifragmentation studies 
[6ll l62l |63| . Due to the small size of our light system, 
this appears to occur for much smaller fragments, and 
even to be reflected in the A = 4 yields. 

The QP formed in Ar + Al, Ar + Ti and Ar + Ni 
collisions are almost identical indicating that we have a 
clean technique for identifying the QP. 
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FIG. 9: (Color online) Average multiplicity of p (a), d (b), t (c), 
3 He (d), a (e) and Li (f) from the QP systems formed in 40 Ar + 
27 Al (open circles), 40 Ti (open triangles) and 58 jVi (solid squares) 
as a function of E* /A. 



V. EXPERIMENTAL EVIDENCE OF CRITICAL 
BEHAVIOR 

We have used several techniques to look for evidence 
of possible critical behavior in the A ~ 36 system. These 
include a Fisher droplet model analysis of the charge dis- 
tributions, searches for the region of maximal fluctua- 
tions using many different observables and tests of the 
fragment topological structure. 



Fisher Droplet Model Analysis of Charge 
Distributions 



The Fisher droplet model has been extensively applied 
to the analysis of multifragmentation since the pioneer- 
ing experiments on high energy proton- nucleus collisions 
by the Purdue group |64l I651 66] . Relative yields of frag- 
ments with 3 < Z < 14 could be well described by a 
power law dependence A~ T suggesting the disassembly 
of a system whose excitation energy was comparable to 
its total binding energy [g^]. The extracted value of the 
power law exponent was 2 < t < 3, which is in a rea- 
sonable range for critical behavior ^j- The success of 
this approach suggested that the multi-fragmentation of 
nuclei might be analogous to a continuous liquid to gas 
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FIG. 10: (Color online) Charge distribution of QP in different 
E* I A window for the reaction 40 Ar + 58 Ni. Lines represent fits. 
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FIG. 11: (Color online) T e f f and X e ff as a function of excitation 
energy for the QP formed in 40 Ar + 58 Ni. 



phase transition observed in more common fluids. 

In the Fisher Droplet Model the fragment mass yield 
distribution, Y(A) , may be represented as 

Y{A) = Y A- T X A2/3 Y A , (10) 

where Yq, t, X and Y are parameters. However, at the 
critical point, X = 1 and Y — 1 and the cluster distri- 
bution is given by a pure power law 

Y(A)=Y A- T , (11) 

The model predicts a critical exponent r ~ 2.21. 

In Fig. we present, for the QP from the reactions 
of 40 Ar + 58 Ni, yield distributions, dN/dZ, observed for 
our nine intervals of excitation energy. 

At low excitation energy a large Z residue always re- 
mains, i.e. the nucleus is basically in the liquid phase 
accompanied by some evaporated light particles. When 
E* I A reaches ~ 6.0 MeV/nucleon, this residue is much 
less prominent. As E*/A continues to increase, the 
charge distributions become steeper, which indicates that 
the system tends to vaporize. To quantitatively pin down 
the possible phase transition point, we use a power law 
fit to the QP charge distribution in the range of Z = 2 - 
7 to extract the effective Fisher-law parameter r e f f by 

dN/dZ - Z^ T "" . (12) 

The upper panel of Fig. 1111 shows r e ff vs excitation en- 
ergy, a minimum with r e // ~ 2.3 is seen to occur in the 
E*/A range of 5 to 6 MeV/nucleon. r e ff ~ 2.3, is close 
to the critical exponent of the liquid gas phase transition 



universality class as predicted by Fisher's Droplet model 
|13|. The observed minimum is rather broad. 

In a lattice gas model investigation of scaling and ap- 
parent critical behavior, Gulminelli et al. have pointed 
out that, in finite systems, the distribution of the maxi- 
mum size cluster, i.e. the liquid, might overlap with the 
gas cluster distribution in such a manner as to mimic the 
critical power law behavior with r e ff ~ 2.2 |67|. They 
further note, however, that at that point the scaling laws 
are satisfied, which suggests a potentially more funda- 
mental reason for the observation of the power law dis- 
tribution |67| . Assuming that the heaviest cluster in each 
event represents the liquid phase, we have attempted to 
isolate the gas phase by event-by-event removal of the 
heaviest cluster from the charge distributions. We find 
that the resultant distributions are better described as 
exponential as seen in Fig. IT21 

The fitting parameter X e ff of this exponential form 
exp{—X e ffZ') was derived and is plotted against excita- 
tion energy in the lower panel of Fig. 1111 A minimum is 
seen in the same region where T e f / shows a minimum. 



B. Maximal fluctuations 

1. Campi Plots 

One of the well known characteristics of the systems 
undergoing a continuous phase transition is the occur- 
rence of the largest fluctuations. These large fluctuations 
in cluster size and density of the system arise because of 
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FIG. 12: (Color online) Same as Fig. llUl but the heaviest cluster 
is excluded on the event by event basis. 



FIG. 13: The Campi plot for different excitation energy windows 
for the QP formed in 40 Ar + 58 Ni. 



the disappearance of the latent heat at the critical point. 
In macroscopic systems such behavior gives rise to the 
phenomenon of critical opalescence |6^| . 

Campi suggested the use of event by event scatter plots 
of the natural log of the size of the largest cluster, lnA max 
versus the natural log of the normalized second moment, 
InSa, of the cluster distribution with the heaviest frag- 
ment removed. For our analysis we use Zmax as the 
measure of the size of the largest cluster and 



S 2 



E 



rii(Zi) 



Zi ■ ni(Zi) 



(13) 



where Zi is the charge number of QP clusters and n.;(Z;) 
is the multiplicity of the cluster Zi. Campi plots have 
proved to be very instructive in previous searches for crit- 
ical behavior jjj. 

In Fig. E| we present such plots for the nine selected 
excitation energy bins. In the low excitation energy 
bins of E*/A < 3.7 MeV/nucleon, the upper (liquid 
phase) branch is strongly dominant while at E*/A > 
7.5 MeV/nucleon, the lower Z max (gas phase) branch is 
strongly dominant. In the region of intermediate E* /A 
of 4.6- 6.5 MeV/nucleon, the transition from the liquid 
dominated branch to the vapor branch occurs, indicating 
that the region of maximal fluctuations is to be found in 
that range. 

Using the general definition of the fcth moment as 



M k = 



iii(Zi). 



(14) 



Campi also suggested that the quantity, 72, defined as 



72 



M 2 M 



(15) 



where Mo, Mi and M2 are the zeroth moment, first 
moment and second moment of the charge distribution, 
could be employed to search for the critical region. In 
such an analysis, the position of the maximum 72 value 
is expected to define the critical point, i.e., the critical 
excitation energy E*, at which the fluctuations in frag- 
ment sizes are the largest. 

The excitation energy dependence of the average values 
of 72 obtained in an event-by-event analysis of our data 
are shown in Fig. 1141 72 reaches its maximum in the 5-6 
MeV excitation energy range. In contrast to observations 
for heavier systems |17|. there is no well defined peak in 
72 for our very light system and 72 is relatively constant 
at higher excitation energies. We note also that the peak 
72 value is lower than 2 which is the expected smallest 
value for critical behavior in large systems. However, 3D 
percolation studies indicate that finite size effects can 
lead to a decrease of 72 with system size [69J, |70j. For 
a percolation system with 64 sites, peaks in 72 under 
two are observed. Therefore, the lone criterion 72 > 2 is 
not sufficient to discriminate whether or not the critical 
point is reached. To carry out further quantitative ex- 
plorations of maximal fluctuations we have investigated 
several other proposed observables expected to be related 
to fluctuations and to signal critical behavior. These are 
discussed below. 
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FIG. 14: (Color online) 72 of the QP systems formed in Ar + 
Al (open circles), Ti (open triangles) and Ni (solid squares) as a 
function of excitation energy. 



2. Fluctuations in the distribution of Z max 

It is supposed that the cluster size distributions should 
manifest the maximum fluctuations around the critical 
point where the correlation length diverges. As a result 
of constraints placed by mass conservation, the size of the 
larg est cluster should then also show large fluctuations 
|68|. Thus, it has been suggested that a possible signal 
of critical behavior is the fluctuation in the size of the 
maximum fragment Recently, Dorso et al. employed 
a molecular dynamics model to investigate fluctuations 
in the atomic number of the heaviest fragment {Z mn .^] 
by determining its normalized variance (NVZ) [2l| , 



NVZ 



(16) 



In that work, they performed calculations of the NVZ 
on two simple systems, one of which should not exhibit 
critical behavior and one which does. For the first they 
used a random partition model in which the population 
of the different mass numbers is obtained by randomly 
choosing values of A following a previously prescribed 
mass distribution [7]]]. In this case the fluctuations in 
the populations arc of statistical origin or are related to 
the fact that the total mass A to t is fixed. No signal of 
criticality is to be expected. In the second case they 
explored the disassembly of systems of the same size em- 
ploying a finite lattice bond percolation model. Such a 
case is known to display true critical behavior [21| . They 
found that that NVZ peaks close to the critical point 
in the percolation model calculation but shows no such 
peak in the random partition model calculation. This 
indicates that the mass conservation criterion, by itself, 
can not induce the peak of NVZ. The details can be 
found in [2l|. 

For our data we plot the normalized variance of 
Zmax/ Zqp as a function of excitation energy in Fig. 1161 
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FIG. 15: (Color online) Z max distributions of the QP systems 
formed in Ar + Al (open circles), Ti (open triangles) and Ni (solid 
squares) in different excitation energy windows. 



A clear maximum, characterizing the largest fluctuation 
of this order parameter, is located in the E* j A ~ 5-6 
MeV/nucleon, 



3. Fluctuations in the distribution of total kinetic energy 

The system which we have studied is a hot system. 
If critical behavior occurs, it should also be reflected in 
large thermal fluctuations. Using a definition similar to 
that of the normalized variance of Z max , we can define 
the normalized variance of total kinetic energy per nu- 
cleon, 



NVE 



(17) 



where E^/A is the total kinetic energy per nucleon and 
(t b m 1 a is its width. Fig. 1171 shows the NVE as a func- 
tion of excitation energy. The observed behavior is very 
similar to that of NVZ. Again, the maximal fluctua- 
tion was found at E* /A = 5 - 6 MeV/u. The maximal 
thermal fluctuations are found in the same region as the 
maximal fluctuations in cluster sizes. 

The use of kinetic energy fluctuations as a tool to mea- 
sure microcanonical heat capacities has also been pro- 
posed [zHzllzl. 
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FIG. 16: (Color online) NVZ of the QP systems formed in Ar + 
Al (open circles), Ti (open traingles) and Ni (solid squares) as a 
function of excitation energy. Vertical line is at 5.6 MeV/u. See 
text. 
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FIG. 17: (Color online) NVE of the QP systems formed in Ar + 
Al (open circles), Ti (open traingles) and Ni (solid squares) as a 
function of excitation energy. Vertical line is at 5.6 MeV/u. See 
text. 



Based on the relation of the heat capacity per nucleon 
and kinetic energy fluctuations, i.e. 
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FIG. 18: (Color online) The kinetical energy fluctuation Aqp jj" 

for the QP formed in 40 Ar + 27 Al (open circles), 48 Ti (open tri- 
angles) and 58 7Vi (solid squares). See text. 



inverting the kinetic equation of state [7£ 



M 



1 \ >^ + 4(M 



i)>r„ 



(19) 



where () indicates the average on the events with the 
same Ej^ and a,i is the level energy density parame- 
ter for fragment i and M is total multiplicity of QP 
particles. A negative heat capacity is indicated if the 
kinetic energy fluctuations exceed the canonical expecta- 
tion A QP (rl/T? n = ck- In 35 MeV/nucleon Au + Au col- 
lisions, the de- excitation properties of an Au QP formed 
at excitation energies from 1 to 8 MeV/nucleon were in- 
vestigated. Abnormal kinetic energy fluctuations were 
observed near the excitation energy previously identified 
as the critical energy [j4| and derived negative heat ca- 
pacities have been taken as a possible signal of the liquid 
gas phase transition. 

In Fig. ^| we present the variable Aqp^^- as a func- 
tion of excitation energy as observed for the present sys- 
tem. The broad peak located at E*/A = 4.0 - 6.5 MeV 
as Fig. El indicates the region of the largest kinetic fluc- 
tuations. We note that the value of this quantity never 
reaches 3/2, which is the canonical expectation. It is 
possible that for such a very small system the finite size 
effects will limit this parameter to values well below the 
canonical expectation as it does for 72. Hence the quanti- 
tative value of heat capacity cy will be difficult to derive 
using Eq. 1181 without a priori knowledge of ck ■ It is 
clear that any value of ck below 0.29 would lead to an 
apparent negative heat capacity. 



where ck and c/ are the kinetic and interaction micro- 
canonical heat capacities per particle calculated for the 
most probable energy partition characterized by a micro- 
canonical temperature T m |7^ . T m can be estimated by 



4- Universal Scaling Laws: A-scalmg 

The recently developed theory of universal scaling laws 
for order-parameter fluctuations has been advanced as 
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providing a method to select order parameters and char- 
acterize critical and non-critical behavior, without any 
assumption of equilibrium |26J. In this framework, uni- 
versal A-scaling laws of the normalized probability dis- 
tribution P[m] of the order parameter to for different 
"system size" (to), should be observed: 



(m) A P[m] = $(Z (A) ) = $[ 



(to) z 



(20) 



with < A < 1, where (to) and to* are the average and 
the most probable values of to, respectively, and <f>(ztA)) 
is the (positive) defined scaling function which depends 
only on a single scaled variable ZtM . If the scaling frame- 
work holds, the scaling relation is valid independent of 
any phenomenological reasons for changing (to) j7|J . The 
A-scaling analysis is very robust and can be studied even 
in small systems if the probability distributions P[m] are 
known with a sufficient precision. 

Botet et al. applied this universal scaling method to 
INDRA data for 136 Xe + 124 Sn collisions in the range 
of bombarding energies between 25 MeV/nucleon and 
50 MeV/nucleon. As the relevant order parameter they 
chose the largest fragment charge, Z max . It was found 
that, at Ei a b > 39 MeV/nucleon, there is a transition in 
the fluctuation regime of Z max . This transition is com- 
patible with a transition from the ordered phase (A = 
1/2) to the disordered phase (A = 1) of excited nuclear 
matter [23j. From this study, they attributed the frag- 
ment production scenario to the family of aggregation 
scenarios which includes both equilibrium models, such 
as the Fisher droplet model, the Ising model, or the per- 
colation model and non-equilibrium models, such as the 
Smoluchowski model of gels. For such scenarios the av- 
erage size of the largest cluster, (Z max ), is the order pa- 
rameter and the cluster size distribution at the critical 
point obeys a power law with r > 2. 

The upper panel in Fig. ^3 shows that A-scaling of 
P[Z m ax) distributions for all E*/A windows above 2.0 
MeV with an assumed A = 1. For our light system, our 
results show that the higher energy data are very well 
scaled with A = 1 (even though not perfectly in the lower 
Z(j±\ tail) but the lower energy data are not. Similar be- 
haviors are also observed in other quantities, such as the 
total kinetic energy per nucleon E^/A fFig. HW bl') and 
the normalized second moment S2 (Fig. EH C )) 01 QP- 
This indicates a transition to A = 1 scaling in the region 
of E*/A = 5.6 MeV. This corresponds to the fluctuations 
of the Z max growing with the mean value, i.e. ^ ma ^ ~ 
constant (see Fig. 12011 . The saturation of the reduced 
fluctuations of Z max (i.e. j§ m ^) observed above corre- 
sponds to the transition to the regime of maximal fluctu- 
ations [24|]. However the lower energy data are not well 
scaled by A = 1/2. 

The pattern of charged fragment multiplicity distribu- 
tions P[n] does not show any significant evolution with 
the excitation energy (Fig. I21f> . and the data are per- 
fectly compressible in the scaling variables of the A = 
1/2 scaling, i.e., the multiplicity fluctuations are small in 
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FIG. 19: (Color online) A-scaling for different quantities: the 
charge distribution of the largest fragment (a), the total kinetic 
energy distribution per nucleon Ej°^/A (b) and the normalized 
second moment (S2) in different E* J A windows. The A=l scaling 
is generally satisfied above 5.6 MeV/nucleon even it is not perfect 
in the lower 2(a) tail. 



the whole excitation energy range. The scaling features 
of experimental Z max ( Fig. EH ( a )) and multiplicity prob- 
ability distributions ( Fig. 1211 ) arc complementary and 
allow one to affirm that the fragment production in Fermi 
energy domain follows the aggregation scenario, such as 
the Fisher droplet model, and two phases of excited nu- 
clear matter with distinctly different patterns of Z max 
fluctuation. It appears that Z max is a very good order 
parameter to explore the phase change |77j . 

From the studies of this section, we conclude that the 
largest fluctuation phase (A = 1) is actually reached 
above 5.6 MeV/nucleon of excitation energy. 



C. Fragment Topological Structure 

In addition to the thermodynamic and fluctuation fea- 
tures of the system, observables revealing some particular 
topological structure may also reflect the critical behav- 
ior for a finite system. For example, if we make a plot 
for the average value of Zi max vs Z max in the different 
excitation energy windows, we immediately see that a 
transition occurs near 5.6 MeV/nucleon (Fig. l22fl . Below 
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FIG. 21: (Color online) A = 1/2 scaling of charged fragment 
multiplicity distributions in different E* /A windows for the QP 
formed in 40 Ar + 58 Ni. 



that point (^moi) increases with decreasing (Z max ). In 
these energy zones, the fragmentation is basically domi- 
nated by evaporation and sequential decay is important. 
But above 5.6 MeV/nucleon excitation energy, (Z^max) 
decreases with decreasing (Z max ). In this region of exci- 
tation, the nucleus is essentially fully vaporized and each 
cluster shows a similar behavior. 

Below we present an exploration of two more detailed 
observables characterizing the topological structure, i.e., 
Zipf law relationships and bimodality. 
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1. Zipf plots and Zipf's law 

Recently, Ma proposed measurements of the fragment 
hierarchy distribution as a means to search for the liq- 
uid gas phase transition a finite system [tsL Fflj . The 
fragment hierarchy distribution can be defined by the so- 
called Zipf plot, i.e., a plot of the relationship between 
mean sizes of fragments which are rank-ordered in size, 
i.e., largest, second largest, etc. and their rank |7ct l79| . 
Originally the Zipf plot was used to analyze the hierar- 
chy of usage of words in a language [$£}, i-e. the relative 
population of words ranging from the word used most 
frequently to the word used least frequently. The integer 
rank was defined starting from 1 for the most probable 
word and continuing to the least probable word. Sur- 
prisingly, a linear relationship between the frequency and 
the order of words was found. Later, many more appli- 
cations of this relationship were made in a broad vari- 
ety of areas, such as population distributions, sand-pile 
avalanches, the size distribution of cities, the distribution 
in strengths of earthquakes, the genetic sequence and the 
market distribution of sizes of firms, etc. It has been sug- 
gested that the existence of very similar linear hierarchy 
distributions in these very different fields indicates that 
Zipf's law is a reflection of self-organized criticality [8l] . 

The significance of the 5-6 MeV region in our data is 
further indicated by a Zipf's law analysis such as that 
proposed in |zl,[7^|. In such an analysis, the cluster size 
is employed as the variable to make a Zipf-type plot, and 
the resultant distributions are fitted with a power law , 



oc rank 



(21) 



where £ is the Zipf's law parameter. In Fig.^Uwe present 
Zipf plots for rank ordered average Z in the nine differ- 
ent energy bins. Lines in the figure are fits to the power 
law expression of Eq. I|21|l . Fig. [21 shows the fitted £ pa- 
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FIG. 23: (Color online) Zipf plots in nine different excitation en- 
ergy bins for the QP formed in 40 Ar + 58 Ni. The dots are data 
and the lines are power-law fits (Eq. 1211 . The statistical error is 
smaller than the size of the circles. 



FIG. 24: (Color online) Zipf parameter as a function of excitation 
energy for the QP formed in 40 Ar + 58 Ni. The position of cross 
illustrates the Zipf law is reached around 5.6 MeV/u excitation 
energy. 



rameter as a function of excitation energy. As shown in 
Fig. 1231 this rank ordering of the probability observa- 
tion of fragments of a given atomic number, from largest 
to the smallest, does indeed lead to a Zipf's power law 
parameter £ = 1 in the 5-6 MeV/nucleon range. When 
£ ~ 1, Zipf's law is satisfied. In this case, the mean size 
of the second largest fragment is 1/2 of that of the the 
largest fragment; That of the third largest fragment is 
1/3 of the largest fragment, etc. 

We note that the nuclear Zipf-type plot which was pro- 
posed in Ref. [T^ nas been applied in the analysis of 
CERN emulsion or Plastic data of Pb + Pb or Plastic 
at 158 GeV/nucleon and it was found that the nuclear 
Zipf law is satisfied when the liquid gas phase transition 
occurs |82^ . 

In a related observation which is consistent with the 
formulation of Zipf's law, percolation model calculations 
[13 suggest that the ratio S p — ^ ma ^ reaches 0.5 
around the phase separation point. Here Zimax is the 
atomic number of the second heaviest fragment in each 
event. Fig. shows S p versus E*/A. S p = 0.5 at 
5.2 MeV/nucleon. It exhibits essentially linear behavior 
(with two different slopes) above and below that point. 

2. Bimodality 

Another proposed test of phase separation is bimodal- 
ity which was suggested in |84||. As has been noted [8{| 
this approach generalizes definitions based on curvature 




<E7A> (MeV) 



FIG. 25: (Color online) The phase separation parameter as a func- 
tion of excitation energy for the QP formed in 40 Ar + 58 Ni. 



anomalies of any thermodynamic potential as a function 
of an observable which can then be seen as an order pa- 
rameter. It interprets a bimodality of the event distribu- 
tion as coexistence, each component representing a differ- 
ent phase. It provides a definition of an order parameter 
as the best variable to separate the two maxima of the 
distribution. In this framework when a nuclear system is 
in the coexistence region, the probability distribution of 
the order parameter is bimodal. 

In analysis of INDRA data jsjj, 
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FIG. 26: (Color online) The average value of bimodal as a function 
of excitation energy for the QP formed in 40 Ar + 58 Ni. 
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-) was chosen as a sorting pa- 



( 

rameter. This parameter may be connected with the 
density difference of the two phases (pi — p g ), which is 
the order parameter for the liquid gas phase transition. 

For our very light system, if we consider the clusters 
with Z < 3 as a gas and the clusters with Z > 4 as a 
liquid, a parameter characterizing the bimodal nature of 
the distribution can be defined as 



P 



E 



Z 4 >4 ' 



<3 ' 



(22) 



Fig.l26lshows the mean value of P as a function of E*/A. 
Here again, the slope shows a distinct change at E*/A = 
5-6 MeV where P — 0, i.e. the point of equal distribution 
of Z in the two phases. 



VI. CALORIC CURVE 

It is also interesting to ask how the caloric curve for 
this light system behaves. Several different experimen- 
tal methods have been applied to the determination of 
caloric curves for nuclear systems. The most common of 
these are the use of slope parameters of the kinetic en- 
ergy spectra or the use of isotopic yield ratios j8|J . Since 
sequential decays and side feeding may be important in 
either case, corrections for such effects must normally be 
applied to observed or "apparent" temperatures in or- 
der to obtain the initial temperatures corresponding to 
the initial excitation energies of the nuclei under investi- 
gation [9(], 0- For heavier systems, a number of mea- 
surements of caloric curves have been reported and 
references therein. In those measurements a flattening 
or plateauing is generally observed at higher excitation 



e4nergies. For light systems such as the A ~ 36 sys- 
tem studied here there are relatively few measurements 
of caloric curves. For construction of the caloric curve 
from the present data we have used both the slope mea- 
surement and isotope ratio technique to derive "initial 
temperatures" from the observed apparent temperatures, 
limiting the use of each to its own range of applicability 
as discussed below. 



A. Low Excitation - the Liquid-Dominated Region 

Determinations based on spectral slope parameters be- 
gan with fitting the kinetic energy spectra for different 
LCPs associated with the nine different bins in excitation 
energy to obtain the apparent slope temperatures T s in 
the QP source frame. T s can be obtained assuming a 
surface emission type Maxwellian distribution, i.e., 
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(23) 
op. 

where E^ in is the kinetic energy in the QP frame and 
V cou i is the barrier parameter. For an example, Fig. I27I 
shows the fits to the kinetic spectra of deuterons and tri- 
tons in the QP frame in four different E* /A windows. 
The dashed lines represent the fits and they show ex- 
cellent agreement with the data. Using such fitted re- 
sults, the excitation function of T s can be obtained for 
each light charged particle as shown in Fig. |2HI In this 
figure, we also plot (dotted line) T s = y / 8E*/A which 
corresponds to the temperature from a simple Fermi gas 
assumption. For these different particles, the apparent 
temperatures are different from each other since the ef- 
fects of sequential decay are different for different par- 
ticles. We note that the temperatures, T s , of 3 He and 
Li, are larger than those of other LCPs, indicating that 
they might be the least affected by the sequential decay 
effects, while T s for protons shows dramatically smaller 
values than the others indicating that the p spectra are 
strongly influenced by later stage emissions. We then 
employed the measured excitation energy dependence of 
the multiplicity for the ejectile under consideration to 
derive initial temperatures from the apparent slope tem- 
peratures mm. 

For each LCP, the measured multiplicity is the sum 
over the entire de-excitation cascade. Since the temper- 
ature of an evaporation residue in an excitation energy 
bin characterized by a small change of excitation from 
Ei to E£ is, to a good approximation, 



(Eini ) 



(M 2 )(T 2 )- (Afi)(Tx 
(M 2 ) - (Mi) 



(24) 



where M 2 and Mi are the multiplicities of a certain LCP 
at the excitation energy and E\ where E\ > E%. The 
details of this method can be found in references |87ll8q |. 

With this method, we can derive the initial tempera- 
tures for each particle. For each particle except protons 
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FIG. 27: (Color online) Kinetic energy spectra in the QP frame 
of 40 Ar + 58 Ni in four selected E* /A windows. Left panels are for 
deuteron and right for tritons. The dots are experimental data and 
the lines are fits with Eq. 1231 
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FIG. 29: (Color online) The deduced caloric curves for the QP 
formed in 40 Ar + 58 Ni. The symbols are displayed in insert. See 
details in text. 



we obtained a reasonable agreement of the respective ini- 
tial temperatures and therefore use their average values, 
as shown by solid squares in Fig. [2HI as a mean initial 
temperature for plotting the caloric curve. For protons, 
the apparent temperature is very low from fits as shown 
in Fig. I28I since a large portion of protons may originate 
from the side feeding besides the sequential decay chain. 
However, the former can not be corrected with Eq. I|24|l . 
It must be emphasized that this technique is based on 
the assumption of sequential evaporation of the ejectiles 
from a cooling compound nucleus source [871 188| . Given 
that the various observables discussed above suggest an 
important transition at 5.6 MeV/u excitation energy, this 
method should not be appropriate above that energy. In 
fact, initial temperatures deduced using this approach ex- 
hibit a very rapid increase at excitation energies above 6 
MeV/nucleon (not shown). We take this as evidence that 
sequential evaporation from a larger parent can not ex- 
plain the multiplicities in the higher energy region. This 
is already suggested by the energy dependence of vari- 
ous multiplicities in Figure 8 as well as by much of the 
discussion in the previous section. Thus we do not em- 
ploy this method based on slope measurements above 5.6 
MeV/nucleon . 



High Excitation (I) - the Vapor-Dominated 
Region 



FIG. 28: (Color online) The apparent temperatures from the 
slopes of different particles (open symbols) and from isotopic ra- 
tio (solid squares) as a function of excitation energy for the QP 
formed in 40 Ar + 58 Ni. The line is the Fermi gas model calcula- 
tion: T s = ^/8E*/A. 



The double isotope ratio temperature technique pro- 
posed by Albergo et al. [8|| has been extensively dis- 
cussed and used in many experiments and theoretical 
calculations. Application of this technique assumes that 
thermal equilibration and chemical equilibration have 
been attained. In an experimental determination, one of 
the major problems is that secondary decay effects can 
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modify the initial temperature strongly 

HUGS EH- 

The experimental apparent double isotope tempera- 
ture, Th~Hb, can be deduced from the ratio of jf^jjt~'- 



14.3 



H—He 



ln[l.m{M d - M a )/(M t - M, He )Y 



(25) 



where M d , M t , Mi. Hei and M a are the isotopic yields 
of d,t, 3 He and a from QP (see Fig.[5J), respectively. As 
shown in Fig. (solid squares), the apparent isotopic 
temperatures are well below those of the simple Fermi gas 
assumption (dotted line in the figure) indicating a strong 
influence due to secondary decay. To estimate the sec- 
ondary decay effect Quantum Statistical Model (QSM) 
calculations were performed to correct the observed dou- 
ble isotope H-He ratio temperatures (Th-Hc) for these 
effects. 

For this purpose we compared results of two different 
calculations, the first published in reference [SH l9lj and 
the second carried out for this work employing the QSM 
model described in reference [8^, El ■ The results of 
the two QSM models are in quite good agreement with 
each other. For nuclei with A ~ 36 in the excitation en- 
ergy range of Interest, averaging results of these models 
indicates that T inlt = (1.75 ± 0.06) x T HHe . The cor- 
rected isotopic temperature is shown in Fig. 1291 as the 
solid circles. 

As emphasized above, this method is based on a model 
which assumes simultaneous fragmentation of a reduced 
density equilibrated nucleus and subsequent secondary 
evaporation from the primary fragments |89|. This 
method should be inappropriate in the lower excitation 
energy where the vapor assumption of the QSM is vio- 
lated. In this case, we do not apply the technique below 
5 MeV/nucleon 

The two techniques differ somewhat in the excitation 
energy range near the transition point, indicating some 
systematic error due to using different techniques in the 
transition point region. This supports the argument for 
restricting the use of each technique to the " appropriate" 
excitation energy region. 

We note that the caloric curve, defined in this manner, 
exhibits no obvious plateau. A polynomial fit to the data 
points leads to a temperature at the transition point of 
8.3 ± 0.5 MeV. 



C. High Excitation (II) - the Ideal Vapor 
Assumption 

If the vapor phase may be characterized as an ideal 
gas of clusters [l3||. then, at and above T = 8.3 MeV, 
this should be signaled by a kinetic temperature, Tki n = 
|E^ n , where E{J n is the Coulomb corrected average ki- 
netic energy of primary fragments. Secondary decay ef- 
fects make it difficult to test this expectation. However, 
in an inspection of the average kinetic energies or appar- 
ent slope parameters (Fig. |^H1 ) for the different species 



observed , we find that, for each excitation energy win- 
dow, the average kinetic energy of 3 He isotropically emit- 
ted in the projectile like frame, is higher than those of 
other species. This together with simple model estimates 
indicates that the 3 He spectra are the least affected by 
secondary decay. Kinetic temperatures for 3 He, defined 
as §(-Efc — B c ) where E)~ is the average kinetic energy and 
B c is the Coulomb energy ( obtained from the fits), are 
plotted as open squares in Fig.|2Hl Above T = 8.3 MeV 
the kinetic temperatures show a similar trend to that of 
the the chemical temperatures but are approximately 1.5 
MeV lower. While not perfect this approximate agree- 
ment provides additional evidence for disassembly of an 
equilibrated system. 

For heavier systems a plateau or flattening is often ob- 
served in caloric curves |l( and the region of entry into 
the plateau appears to be very close to the point which 
has been identified as the point of maximal fluctuations. 
The reason for this flattening is still under discussion. 
It may reflect expansion and/or s pin odal decom posi tion 
inside the coexistence region pi fall l9l I9I I9I l96j . In 
contrast, our light system does not show the flattening. 
This suggests which that the transition under investiga- 
tion may differ from that seen in the heavier systems. 

Taken together with the observations indicating max- 
imal fluctuations and the particular features of the frag- 
ment topological structure at 5.6 MeV/u excitation, the 
comportment of this caloric curve provides further evi- 
dence suggesting that the observed transition is taking 
place at, or very close to, the critical point. 



VII. DETERMINATION OF CRITICAL 
EXPONENTS 

Since the pioneering work on extraction of the critical 
exponents for nuclear multifragmention from EOS data 
[93, several additional experimental and theoretical ef- 
forts have been attempted [H |H [l0(| . In the latter 
works, Elliott et al. show that the scaling behavior can 
remain even in small systems and the critical exponents 
can be extracted. 

In the Fisher droplet model, the critical exponent r 
can be deduced from the cluster distribution near the 
phase transition point. In Sec.V(A), we already deter- 
mined, from the yield distributions, r e // ~ 2.31 ± 0.03, 
which is close to that for the liquid gas phase transition 
universality class. In terms of the scaling theory, t can 
also be deduced from, (S corr ), the slope of the correla- 
tion between In (S3) vs ln(S 2 ) (HI, where S3 = M3/M1, 
shown in Fig. [301 is related to t as 



35, 



1 



(26) 



Assuming the value of T c — 8.3 MeV as determined from 
our caloric curve measurements, we explored the correla- 
tion of of S2 and S3 in two ranges of excitation energy C 
see Figure 31. The moments were computed by exclusion 
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FIG. 30: (Color online) The correlation between ln(S3) vs ln(S2) 
and the linear fit. 



FIG. 31: (Color online) 

Zmax as a function of T/T c (a) and the 
extraction of the critical exponent /3 (b). 



of the species with Z max in the " liquid" phase but inclu- 
sion in the "vapor" phase. The slopes were determined 
from linear fits to the "vapor" and "liquid" regions re- 
spectively and then averaged. In this way, we obtained a 
value of t= 2.13 ±0.1. See Fig. 

Other critical exponents can also be related to other 
moments of cluster distribution, which were defined 
in Eg. 1141 Since, for our system, we have already deduced 
the initial temperatures and determined a critical tem- 
perature T c — 8.3 MeV at point of maximal fluctuations, 
we can use temperature as a control parameter for such 
determinations. In this context, the critical exponent (3 
can be extracted from the relation 

Z max oc (1 - (27) 

1 c 

and the critical exponent 7 can be extracted from the 
second moment via 

M 2 oc|l-^r. (28) 

In each, |1 — ^-\ is the parameter which measures the 
distance from the critical point. 

Fie. BT1 explores the dependence of Z max on (ijr). We 
note a dramatic change of Z max around the critical tem- 
perature T c . LGM calculations also predict that the slope 



of Z m nx vs T will change at the liquid gas phase tran- 
sition 101] . Physically, the largest fragment is simply 
related to the order parameter pi — p g (the difference 
of density in nuclear 'liquid' and 'gas' phases). In infi- 
nite matter, the infinite cluster exists only on the 'liquid' 
side of the critical point. In finite matter, the largest 
cluster is present on both sides of the phase transition 
point. In this figure, the significant change of the slope 
of Z max with temperature should correspond to a sudden 
disappearance of the infinite cluster ('bulk liquid') near 
the phase transition temperature. For the finite system, 
it reflects the onset of critical behavior there. Using the 
left side of this curve (i.e. liquid side), we can deduce the 
critical exponent (3 by the transformation of the x axis 
variable to the distance from the critical point. Fig. 131b 
shows the extraction of (3 using Eq. |23 An excellent fit 
was obtained in the region away from the critical point, 
which indicates a critical exponent (3 = 0.33 ± 0.01. Near 
the critical point, the finite size effects become stronger 
so that the scaling law is violated. The extracted value 
of (3 is that expected for a liquid gas transition (See Ta- 
ble.II) HI. 

To extract the critical exponent 7, we take M2 on the 
liquid side without Z max but take Mi on the vapor side 
with Z max included. Fig. 1321 shows ln(M2) as a function 
of ln(|l — tjt\). The lower set of points is from the liq- 
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FIG. 32: (Color online) The extraction of critical exponents 7. 
See texts for details. 



uid phase and the upper set of points is from the vapor 
phase. For the liquid component, we center our fit to 
Eq.EHlabout the center of the range of (1 — T/T c ) which 
leads to the linear fit and extraction of (3 as represented 
in Figure. I3*T1 We obtain the critical exponent 7 = 1.15 
± 0.06. This value of 7 is also close to the value expected 
for the liquid gas universality class (see Table II). It is 
seen that the selected region has a good power law de- 
pendence. However, a similar effort to extract the 7 in 
the gas phase is not successful: a small value less than 

0. 20 is deduced. This may be due to the finite size ef- 
fects for this very light system. Since the largest cluster 
still exists in the vapor side, its inclusion (or exclusion) 
in M2 might perturb the determination of the moment, 
resulting in an imprecise value of 7 extracted from the 
vapor phase. For comparison, we just show, for the vapor 
phase, a line representing the 7 derived from the liquid 
side. This line only agrees with the last few vapor points, 

1. e. the highest temperature points (the contamination 
of M 2 should the least there) . 

Since we have the critical exponent (3 and 7, we can 
use the scaling relation 



1 



£ + 7 



(29) 



to derive the critical exponent a. In such way, we get the 
(7 = 0.68 ± 0.04, which is also very close to the expected 
critical exponent of a liquid gas system. 

Finally, it is possible to use the scaling relation 



(3 



(30) 



to check the r value which was determined from the 
charge distributions using Fisher droplet model power 
law fits around the critical point (see Fig. I10|) . Using 
Eq. EH we Obtain r = 2.22 ± 0.46, which, though less 
precise, is in agreement with the values of 2.31 ± 0.03 ob- 
tained from the charge distribution around the point of 



TABLE II: Comparison of the Critical Exponents 
Exponents 3D Percolation Liquid-Gas This work 



7 
a 



2.18 

0.41 

1.8 

0.45 



2.22±0.46 (Eq.dJ 
2.21 2.31T0.03 (Eg. [HI 

2.13T0.10 (Eq.HU 
0.33 0.33T0.01 
1.23 1.15T0.06 
0.64 0.68±0.04 



maximal fluctuations and 2.15 ± 0.1 extracted from the 
correlation of ln^) vs ln^). 

To summarize this section, we report in Table. II a com- 
parison of our results with the Values expected for the 
3D percolation and liquid gas system universality classes 
and with the results obtained by Elliott et al. for a heav- 
ier system. Obviously, our values for this light system 
with A~36 are consistent with the values of the liquid 
gas phase transition universality class rather the 3D per- 
colation class. 



VIII. CONCLUSIONS 

In conclusion, an extensive survey of the features of 
the disassembly of nuclei with A ~ 36 has been reported. 
To carry out this analysis, the de-excitation products of 
the A~ 36 quasi projectile source were first reconstructed 
using a new technique based upon on the three source fits 
to the light particle spectra and use of a rapidity cut for 
IMF. Monte Carlo sampling techniques were applied to 
assign all particles to one of three sources (QP, NN and 
QT). 

At an excitation energy ~ 5.6 MeV/nucleon key ob- 
servables demonstrate the existence of maximal fluctua- 
tions in the disassembly process. These fluctuation ob- 
servables include the Campi scattering plots and the nor- 
malized variances of the distributions of order param- 
eters, (Z max ) and total kinetic energy. Recently pro- 
posed A-scaling analysis also show a universal behavior 
at higher excitation energy where the saturation of the 
reduced fluctuations of Z max (i.e. j§ m ^ s y) is observed. 
This corresponds to the transition to a regime of large 
fluctuations from an ordered phase at lower excitation 
energy. 

At the same excitation energy ~ 5.6 MeV/nucleon, the 
Fisher droplet model prediction is satisfied, with a Fisher 
power law parameter, r = 2.3, close to the critical expo- 
nent of the liquid gas phase transition universality class. 
In addition, the fragment topological structure shows 
that the rank sorted fragments obey Zipf 's law, proposed 
as a signature of liquid gas phase transition [73, at the 
maximal fluctuation point. The related phase separa- 
tion parameter |8jj] shows a significant change of slope 
with excitation energy. The correlation of the heaviest 
fragment and the second heaviest fragment demonstrates 
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a transition around 5.6 MeV/u of excitation energy. A 
bimodality test [84| also gives an indication of a phase 
change in the same excitation energy region. 

The caloric curve shows a monotonic increase in tem- 
perature and no plateau region is apparent, in contrast to 
caloric curves seen for heavier systems Q . At the appar- 
ent critical excitation energy the temperature is 8.3 ±0.5 
MeV. Taking this to be the critical temperature for this 
system, we extracted the critical exponents (3, 7 and a. 
The deduced values are consistent with the values of the 
liquid gas phase transition universality class 68] . 

Since some fluctuation observables, such as the struc- 
ture of the Campi plot, Z max fluctuations etc., could be 
produced by mass conservation effects where no assump- 
tion of the critical behavior is needed, these observables 
by themselves do not guarantee that the critical point 
has been reached. What differentiates the present work 
from previous identifications of points of critical behav- 
ior in nuclei, in addition to the fact that these are the 
lightest nuclei for which a detailed experimental analysis 
has been made, is the comportment of the caloric curve 
and the critical exponent extraction. Taken together, 



this body of evidence suggests a liquid gas phase change 
in an equilibrated system at, or extremely close to, the 
critical point. Detailed theoretical confrontations with 
models which include or exclude a liquid gas phase tran- 
sition are certainly interestin g and we lcome. Some work 
along this line is in progress |l02l Il03j . 
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